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Abstract 

The act of measuring a physical signal or field suggests a generalization of the wavelet 
transform that turns out to be a windowed version of the Radon transform. A reconstruc- 
tion formula is derived which inverts this transform. A special choice of window yields 
the "Analytic-Signal transform" (AST) , which gives a partially analytic extension of func- 
tions from R n to C n . For n =1, this reduces to Gabor's classical definition of "analytic 
signals." The AST is applied to the wave equation, giving an expansion of solutions in 
terms of wavelets specifically adapted to that equation and parametrized by real space and 
imaginary time coordinates (the Euclidean region). 

1. Introduction 

The ideas presented here originated in relativistic quantum theory [13, 14, 15], where 
a method was developed for extending arbitrary functions from R n to C n in a semi- 
analytic way. This gave rise to the "Analytic-Signal transform" (AST) [16]. Later it was 
realized that the AST has a natural generalization to what we have called a Windowed 
X-Ray transform [17], and the latter is a special case of a Windowed Radon transform, 
to be introduced below. For n = 1, these transforms reduce to the (continuous) Wavelet 
transform. In the general case, they retain many of the properties of the Wavelet transform. 

In Section 2 we motivate and define the d-dimensional Windowed Radon transform in R n 
for 1 < d < n and derive reconstruction formulas which can be used to invert it. In Section 
3 we define the AST in R n and give some of its applications. In Section 4 we develop a 
new application of the AST by generalizing a construction in [16] to the wave equation 
m R 2 . This results in a representation of solutions of the wave equation as combinations 
of "dedicated" wavelets that are especially customized to that equation. In particular, 
these wavelets are themselves solutions and represent coherent wave packets, being well- 
localized in both space (at any particular time) and frequency, within the limitations of the 
uncertainty principle. The parameters labeling these wavelets (i.e., the variables on which 
the AST depends) have a direct geometrical significance: They give the initial position, 
direction of motion and average frequency or color of the wavelets. The representation of a 
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solution in terms of these wavelets therefore gives a geometrical-optics (ray) picture of the 
solution. It is suggested that this could be of considerable practical value in signal analysis, 
since many naturally occurring signals (e.g., sound waves, electromagnetic waves) satisfy 
the wave equation away from sources and the geometrical-optics picture gives a readily 
accessible display of their informational contents. 

2. Windowed Radon Transforms 
2.1. The Windowed X— Ray Transform 

Suppose we wish to measure a physical field distributed in R n . This field could be a 
"signal," such as an electromagnetic field or the pressure distribution due to a sound wave. 
For simplicity, we assume to begin with that it is real-valued, such as pressure. (Our 
considerations easily extend to complex- valued, vector- valued or tensor- valued signals, 
such as electromagnetic fields; we shall indicate later how this is done.) The given field is 
therefore a function / : R n — > R. We may think of R n as physical space (so that n = 3), 
in which case the field is time-independent, or as space-time (so that n = 4), in which case 
the field may be time-dependent. In the former case, R n is endowed with a Euclidean 
metric, while in the latter case the appropriate metric is Lorentzian, as mandated by 
Relativity theory. 

Actual measurements are never instantaneous, nor do they take place at a single point 
in space. A measurement is performed by reading an instrument, and the instrument 
necessarily occupies some region in space and must interact with the field for some time- 
interval before giving a meaningful reading. Let us assume, to begin with, that the spatial 
extension of the instrument is negligible, so that it can be regarded as being concentrated at 
a single point at any time. We allow our instrument to be in an arbitrary state of uniform 
motion, so that its position is given by x(t) = x+vt, where i G Risa "time" parameter and 
x, v G R n . Note that t need not be the physical time. For example, if R n is space-time, 
then each "point" x represents an event, i.e. a particular location in space at a particular 
time. In that case, the line x(£) is called a world-line and represents the entire history of 
the point-instrument. The "velocity" vector v then has one too many components and 
may be regarded as a set of homogeneous coordinates for the physical velocity. Note that 
in this case v cannot vanish, since this would correspond to an instrument not subject to 
the flow of time. Even if R n is space, the case v = is not interesting since then the 
instrument can only measure the field at a single point. We therefore assume that v^O, 
hence v e R* = R n \{0}. 

Let us assume that the reading registered by the instrument at time s gives a weight h(t — s) 
to the value of the field passed by the instrument at time t. (For motivational purposes 
we note that causality would demand that h(t — s) = for t > s; moreover, h(t — s) 
should be concentrated in some interval s — r < t < s, where r is a "response time" or 
memory characteristic of the instrument. However, the results below do not depend on 
these assumptions.) Our model for the observed value of the field at the "point" x, as 
measured by the instrument traveling with uniform velocity v, is then 
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/oo 
dth(t)f(x + vt). 
-oo 



(1) 



To accomodate complex-valued signals, we allow the weight function h to be complex- 
valued. h(t)* will denote the complex conjugate of h(t). In order to minimize analytical 
subtleties, we assume that h is smooth and bounded, and that / is smooth with rapid 
decay (say, a Schwartz test function). 

Definition 1. The Windowed X-Ray Transform of f : R n — > C is tie function f : 
R n x R? -> C given by 

/oo 
dth(t)*f(x + vt). (2) 
-oo 



Remarks. 

1. In the special case h(t) = 1 and |v| = 1, is known as the (ordinary) X-Ray transform of 
/ (Helgason [11]), due to its applications in tomography. We may then regard fh as being 
defined on the set of all lines in R n , independent of their parametrization. In the general 
case, we think of the function h(t) as a window, which explains our terminology. 

2. Some work along related lines was recently done by Holschneider [12]. He considers a two- 
dimensional wavelet transform which is covariant under translations, rotations and and 
dilations of R 2 . When the window function is supported on a line, say h(ti,t2) = #(£2), 
this becomes an X-Ray transform in R 2 . His inversion method is less direct than ours in 
that it involves a limiting process. 

3. Note that fh has the following dilation property for a^O: 

/oo 
dt la)" 1 h(t/a)* /(x + tv) = f ha (x, v) (3) 
-00 

where h a (t) = |a| _1 h(t/a). This may be used to study the behavior of fh as v — > 0. For 
the "forbidden" value v = 0, the transform becomes /^(x, 0) = /i(0)*/(x), where h is the 
Fourier transform of h. (We shall see that h(0) = for "admissible" h.) 

4. For n = 1 and v 7^ 0, a change of variables gives 

= \v\- 1 l 2 Wf(x,v), 

where W f is the usual wavelet transform of / [5, 7, 23], with v playing the role of a dilation 
factor and the window function h(t) playing the role of a basic wavelet. 

5. All our considerations extend to vector-valued signals. The cleanest approach is to let the 
window function h* assume values in the dual vector space, so that h(t)*f(x + vt) and 
//i(x, v) are scalars. More than one window needs to be used (or rotated versions of a 
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single window), in order to 'probe' the different components of /. The same applies to 
tensor-valued signals such as electromagnetic fields, since they may be regarded as being 
valued in a higher-dimensional vector space. However, a more correct way to measure a 
vector- or tensor-valued field is to use an instrument which is not rotationally invariant, 
and that implies that the instrument has some spatial extension. This is done in Section 
2.3. 

It will be useful to write fh in another form by substituting the Fourier representation of 
/ into fh- Formally, this gives 

/ h (x,v)= [°° dt [ dpe 2 ^ +t ^h(tyf(p) 

J -oc iR" 

= I dpe 2 ^ x Mp-v)*/(p) 
</r™ 

= (^x,v,/)z, 2 (dp) = (^x,v,/)i2( dx ), 



(5) 



where /i XjV is defined by 



so that 



h^ v (p)=e- 2 ^h(p-v), (6) 



/i x ,v(x') = / dpe 2 " p -( x '- x ^(p-v). (7) 

(We have adopted the convention used in the physics literature, where complex inner 
products are linear in the second factor and antilinear in the first factor.) 
The functions /i X)V are n-dimensional "wavelets" and will be used in the next subsection 
to reconstruct the signal /. Note that /i XjV (hence also /i X)V ) is not square-integrable for 
n > 1, since its modulus is constant along directions orthogonal to v. But eq. (5) still 
makes sense provided / is sufficiently well-behaved. (This is one of the reasons we have 
assumed that / is a test function.) 

A common method for the construction of n-dimensional wavelets consists of taking tensor 
products of one-dimensional wavelets. However, this means that not all directions in R n 
are treated equally, and consequently the set of wavelets does not transform "naturally" 
(in a sense to be explained below) under the affme group G of R n , which consists of all 
transformations of the form 

x g(A, b)x = Ax + b (8) 

with A a non-singular n x n matrix and b e R n . Each such g(A, b) defines a unitary 
operator on L 2 (R n ), given by 

(^(Aty/Hx^i^ri/^-^x-b)), (9) 

where \A\ denotes the absolute value of the determinant of A. The map 
g(A,b) i— > U(A,h) forms a representation of G on L 2 (R n ), meaning that it preserves 
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the group structure of G under compositions. To see how /i x v transforms under U, note 
that the unitarity of U implies 



(U(A,b) h^,f) = {h^,U(A,b)- 1 f) 

/oo 
dth(t)*\A\$f(A(x + tv) + b) (10) 
-oo 

= \A\*( tlAx+h,Av , f )• 

Hence 

U(A,b)h x ,v=\A\*h Ax+bjA v, (11) 

which states that affine transformations take wavelets to wavelets. Thus, for example, 
translations, rotations and dilations merely translate, rotate and dilate the labels {x, v}, 
while the factor \A\? preserves unitarity. By contrast, tensor products of one-dimensional 
wavelets are not transformed into one another by rotations. 

2.2. A Reconstruction Formula 

A reconstruction consists of a recovery of / from f h or its restriction to some subset. In 
the one-dimensional case, for example, / can be reconstructed using all of Rx R* or (for 
certain choices of h) just a discrete subset [2, 6, 18, 21, 22]. For general n, the choice 
of reconstructions becomes even richer since various new possibilities arise. For example, 
h may have symmetries which imply that fh is determined by its values on some lower- 
dimensional subsets of R n x R™, making integration over the whole space unnecessary 
and, moreover, undesirable since it may lead to a divergent integral. Furthermore, / may 
satisfy some partial differential equation which implies that it is determined by its values 
on subsets of R n . For example, if R n is space-time and / represents a pressure wave or 
an electromagnetic potential, it satisfies the wave equation away from sources, hence is 
determined by initial data on a Cauchy surface in R n , and it becomes both unnecessary 
and undesirable to use all of R n x R™ in the reconstruction (cf. Sections 3.2 and 4). 
The reconstruction to be developed in this subsection is "generic" in that it does not 
assume any particular forms for h(t) or /(x). It uses all of R n x R™, so it breaks down 
for certain choices of h or /. Again we emphasize that this is far from the only way to 
proceed; other types of reconstruction will be discussed below and elsewhere. The present 
reconstruction formula is interesting in part because it generalizes the one for the ordinary 
continuous wavelet transform (n = 1). 

To reconstruct /, we look for a resolution of unity in terms of the vectors /i x ,v This means 
we need a measure <i//(x, v) on R n x R™ such that 

/ <fc(x,v)|/ fc (x,v)| 2 = / rfx|/(x)| 2 =||/||l 2 . (12) 

(Such an identity is sometimes called a "Plancherel formula.") For then the map T: f i— > fh 
is an isometry from L 2 (dx.) onto its range H. C L 2 (d[i), and polarization gives 

(g,T*Tf) LHdx) = {Tg,Tf) n = (gJ) LHdx) . (13) 
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This shows that / = T*Tf = T* fh in L 2 (<ix), which is the desired reconstruction formula. 
(Cf. [16] for background on resolutions of unity, generalized frames and related subjects.) 
To obtain a resolution of unity, note that 

/ fc (x,v)= (/*(p-v)*/(p))~(x), (14) 
where " denotes the inverse Fourier transform, so by Plancherel's theorem, 

/ rfx|A(x,v)| 2 = / rfp|Mp-v)| 2 |/(p)| 2 . (15) 
Jr.™ Jr» 

We therefore need a measure dp{\) on R™ such that 

H(p) = / dp(v) \h(p • v)| 2 = 1 for almost all p, (16) 

since then d/i(x, v) = dx.dp(v) has the desired property. The solution is simple: Every 
p ^ can be transformed to q = (1, 0, • • • , 0) by a dilation and rotation of R n . That is, 
the orbit of q (in Fourier space) under dilations and rotations is R™. Thus we choose dp 
to be invariant under rotations and dilations, which gives 

dp(y) = N\v\~ n dv, (17) 

where N is a normalization constant, |v| is the Euclidean norm of v and dv is Lebesgue 
measure in R n . Then for p ^ 0, 



H(p) = H(<i) = N J \v\- n dv\h( Vl )\ 2 

Now a straightforward computation gives 

dv 2 ■ ■ ■ dv n 7T n / 2 



(18) 



L 



! Rn -i {vl + ---vl)^ |«i|r(n/2)- (19) 

This shows that the measure dp(x, v) = dx.dp(v) gives a resolution of unity if and only if 

Oh = H0\ 2 <oo, (20) 

which is precisely the admissibility condition for the usual (one-dimensional) wavelet 
transform [5]. (As mentioned above, admissibility implies that h(0) = 0.) If h is admissible, 
the normalization constant is given by 

T(n/2) , s 

N = \' ' 21 

7T n / 2 C h V ' 

and the reconstruction formula is 
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Z(x') = (T*A)(x') = N f |v|-»dxdv/i XiV (x') A(x, v). (22) 

</R™xR™ 

The sense in which this formula holds depends on the behavior of /. The class of possible 
/'s, in turn, depends on the choice of h. Note that in spite of the factor |v| n in the denom- 
inator, there is no problem at v = since /^(x, 0) = h(0)* /(x) = by the admissibility 
condition, and a similar analysis can be made for small |v| by using the dilation property 
(eq.(3)). 

2.3. The d— Dimensional Windowed Radon Transform 

Next, we allow the instrument to extend in k > spatial dimensions. For example, k = 1 
for a wire antenna whereas k = 2 for a dish antenna. If R n is space, then k < n; if 
R n is space-time, then k < n — 1. When moving through space with a uniform velocity, 
the instrument sweeps out a <i-dimensional surface in R n , where d = k + 1 if the motion 
is transverse to its spatial extension and d = k if it is not. If k < n, then the set of 
non-transversal motions is "non-generic" (has measure zero) and can thus be ignored; we 
therefore set d = k + 1 in that case. If k = n, then necessarily d = n. In either case, we 
represent the moving instrument by a window function h : R d — > C. 

The parameter t E R has thus been replaced by t e R d . The velocity vector v, which 
may be regarded as a linear map 1 1— > vt from R to R n , is now replaced by a linear map 
A : K d -> R n , which we call a motion of the instrument in R n . Denote the set of all 
such maps by L(R d , R n ). Later, when seeking reconstruction, we shall need to restrict 
ourselves to subsets of L(R d , R n ) ('rigid' motions); but this need not concern us presently. 

Definition 2. The d-dimensional Windowed Radon Transform of f : R n — > C is the 
function f h :R n x L(R d , R n ) -> C given by 

fh(x,A)= [ dth(t)*f(x + At). (23) 

Upon substituting the Fourier representation of /, a computation similar to the above 
yields the expression 

f h (x,A)= [ dt [ dpe 2 ^ +A ^h(tyf(p) 

= / dpe 27Tip -*h(A' P y /(p) ( 24 ) 

= (hjc,A,f)L 2 (dp) = (h x ,A, f)L 2 {dx), 

where A' : R n — > R d is the map dual to A. (For given bases in R d and R n , A is represented 
by an n x d matrix; then A' is the transposed d x n matrix.) In the above equation we 
have set 

K A {p)=e- 2 ^h(A l p), (25) 
which gives the generalized wavelets 
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fcx,A(x')= / dpe 2 ™»<*'-^h(A'p). 



(26) 



Let us now attempt to reconstruct / from fh by generalizing the procedure in Section 2.2. 
Eq. (15) now becomes 

/ dx|/,(x,A)| 2 = f dp|^(A'p)| 2 |/(p)| 2 . (27) 

Again we need a measure dp(A) which is invariant under dilations and rotations of R n . 
Now the largest set of maps A which can be considered consists of all those with rank d. 
(Otherwise the instrument sweeps out a surface of dimension lower than d.) Let us call 
this set L d (R d ,R n ). Then a measure on L d (R d , R n ) which is invariant with respect to 
rotations and dilations of R n has the form 

d~p{A) = | det(A'A) | " n / 2 dA, (28) 

where dA is the Haar measure on L(R d , R n ) ss R nd as an additive group. However, 
no reconstruction is possible using the measure d^dp(A) on R n x L d (R d ,R n ), because 
no admissible window exists in general when d > 1. (This can be easily verified when 
d = n = 2.) Thus L d (R d , R n ) is too large, and we return to our imaginary measuring 
process for inspiration. On physical grounds, we are interested in rigid motions of the 
instrument. A map corresponding to such a motion must have the form A = vRJ, where 
J : R d — > R n is the canonical inclusion map, R is a rotation of R n (RJ gives the orientation 
of the instrument as well as its direction of motion), and v > is the speed. (If R n is space- 
time, then "rotations" involving the time axis are actually Lorentz transformations! For the 
present, assume that R n is space, so R is a true rotation.) We therefore parametrize the set 
of permissible A : s by (v, R) G R + x SO(n) = G, where SO(n) is the group of unimodular 
orthogonal n x n matrices, which represent rotations in R n . This parametrization is 
redundant because two rotations of R n which have the same effect on the subspace R d 
give the same motion. A non-redundant parametrization of rigid motions is given by 
R + x (SO(n)/SO(n — d)). However, we use the redundant one here for simplicity. (We 
shall need the Haar measure on SO(n).) Note that for d = 1, J is represented by the 
vector q = (1, 0, ■ • ■ , 0) and the set of all maps A = vRJ as above coincides with the set 
R" of non-zero velocities considered in Sections 2.1 and 2.2. A measure on G which is 
invariant under rotations and dilations of R n (i.e., under G itself) has the form 

dp{A) = Nv~ l dvdR, (29) 

where is a normalization constant and dR is the Haar measure on SO(n). Thus for all 
Pt^O, 

ff(p)= / dp(A)\h(A'p)\ 2 = H(ci). (30) 
Jg 

Now 
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A'ci = vJ'R , q = vJ , -R' 1 , (31) 

where Ri is the first row of R, which is a unit vector, and J' is the projection of R n onto 
R d . The admissibility condition therefore reads 

N' 1 = [ v^dv [ dR \h(v J'K'^l 2 <oo. (32) 

J0 JSO(n) 

For d = 1, this reduces to eq. (20). If h is admissible, we obtain the reconstruction formula 

/(x') = / rfx / dp(A)h x , A (x') f h (pc,A). (33) 



3. Analytic— Signal Transforms 
3.1. Analytic Signals in One Dimension 

Suppose we are given a (possibly complex-valued) one-dimensional signal 
/ : R — > C. For simplicity, assume that / is smooth with rapid decay. Consider the 
positive- and negative- frequency parts of /, defined by 



oo 

2tt ipx 



f + (x)= I dpe^ x f(p) 

,o (34) 



f~(x)= [ dpe 2 ^ x f(p), 

J —oo 



Then /+ and / extend analytically to the upper-half and lower-half complex planes, 
respectively; i.e., 



/•OO 

f+(x + iy)= / dpe 27Tip ( x+iy) f(p), y>0 
Jo 

f-(x + iy)= [ dpe 2 ^ x+i ^f(p), y<0, 

J — oo 



(35) 



since the factor e 27Tpy decays rapidly for p — > ±oo in the respective integrals. f + (z) 
and f~(z) are just the (inverse) Fourier-Laplace transforms of the restrictions of / to the 
positive and negative frequencies. If / is complex- valued, then /+ and / _ are independent 
and the original signal can be recovered from them as 

f(x) = lim [f + (x + iy) + f~(x - iy)] . (36) 
yio 

If / is real-valued, then 

f(p) = f(-p)*- (37) 
9 



In that case, /+ and / are related by reflection, 

f + (x + iy) = f~(x - iy)*, y > 0, 



and 



f(x) = 2 lim$tf + (x + iy)=2 lim£/-(a; - iy). 

ylO ylO 



(38) 



(39) 



When / is real, the function f + (z) is known as the analytic signal associated with f(x). 
Such functions were first introduced and applied extensively by Gabor [8]. A complex- 
valued signal would have two independent associated analytic signals / + and /~. What 
significance do f ± have? For one thing, they are regularizations of /. Eq. (36) states that 
/ is jointly a "boundary-value" of the pair / + and / _ . As such, / may actually be quite 
singular while remaining the boundary-value of analytic functions. Also, 2f^ provide a 
kind of "envelope" description of / (cf. Born and Wolf [4], Klauder and Sudarshan [19]). 
For example, if f(x) = cosx, then 2f ± (x) = e ±lx . 

In order to extend the concept of analytic signals to more than one dimension, let us first 
of all unify the definitions of /+ and /~ by setting 



f(x + iy) 



dp9{py)e 2mp{x+iy) f{p) 



(40) 



for arbitrary x + iy G C, where 9 is the unit step function, defined by 





i: 


u < 


6{u) = < 


2' 


u = 




,1, 


u > 



(41) 



Then we have 



f + (x + iy), 
f(x + iy) = <j y(x), 

, f~(x + iy), 



y>0 
y = 
y<0. 



(42) 



Although this unification of / + and / _ may at first appear to be somewhat artificial, it 
turns out to be quite natural, as will now be seen. Note first of all that for any real u, we 
have 



9{u)e 



-2iru 



1 dr 

2ni J_ 00 t - i 



(43) 



since the contour on the right-hand side may be closed in the upper-half plane when u > 
and in the lower-half plane when u < 0. For u = 0, the equation states that 
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i 

2^ 



oo 
oc 



dr 



(44) 



r 



2 + l 2' 



in agreement with our definition, if we interpret the integral as the limit when L — > oo of 
the integral from — L to L. Therefore 



9(py) e 



1 



dr 



3 27rip(x+rj/) 



2m ./_oo r - / 



(45) 



If this is substituted into our expression for f(z) and the order of integrations on r and p 
is exchanged, we obtain 



f ~ {x+w) = 2^I_ 00 7h f{x+Ty) 



(46) 



for arbitrary x + iy G C. We have referred to the right-hand side as the Analytic-Signal 
transform of f(x) [16, 17]. It bears a close relation to the Hilbert transform, which is 
defined by 



1 POO 7 

Hf(x) = -FV / —/(*-«), 



(47) 



where PV denotes the principal value of the integral. Consider the complex combination 

1" 



f( x )- iH f(x) = - I 

J- c 



du 



mS(u) + PV- 



u 



f(x-u) 



i r du 

— lim / j(x — u) 

Jfi ej.0 J _ oc u — it 

— lim / ^ T fix — rt) 
7ri eio /_ m r - % 



(48) 



Similarly, 



Hence 



= 2 lim f (a; — it). 



/(x) + = 2 lim f(x + it). 

ej.0 



Hf(x) = i lim [/(a; - it) - f(x + it)], 

ej.0 



(49) 



(50) 



which for real-valued / reduces to 



Hf(x) = 2 lim $sf(x + it) = —2 lim Qf(x — it). 

ej.0 ej.0 



(51) 
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3.2. Generalization to n Dimensions 

We are now ready to generalize the idea of analytic signals to an arbitrary number of 
dimensions. Again we assume initially that /(x) belongs to the space of Schwartz test 
functions £(R n ), although this assumption proves to be unnecessary. 

Definition 3. The Analytic-Signal Transform (AST) of f e S(R n ) is the function f : 
C n -> C defined by 

W..;=* /(x+Ty) - (52) 

The same argument as above shows that for z = x + iy e C n , 

/(z)= / d n P e(p-y)e 2 ^ z f(p) 

Mn (53) 



= f dpe 2 ^f(p) 

JM y 

where M y is the half-space 



M y = {p G R n such that p • y > 0}, y ^ 0. (54) 

We shall refer to the right-hand side of eq. (53) as the (inverse) Fourier-Laplace transform 
of / in My. The integral converges absolutely whenever / e L 1 (R n ), since |e 27 " p ' z | < 1 on 
M y , defining / as a function on C n , although not an analytic one in general (see below). 
This shows that /(z) can actually be defined for some distributions /, not only for test 
functions. 

Note: In spite of the appearance of expressions such as p • y, we have not assumed any 
particular metric structure in R n . The Fourier transform naturally takes functions on 
R n (considered as an abelian group) to functions on the dual space R n = (R n )* of linear 
functionals, and p-y merely denotes the value p(y). (See Rudin [25].) This remark becomes 
especially important when considering time-dependent signals, so that R n is space-time, 
for then the natural structure on R n is a Lorentzian metric rather than a Euclidean metric 
(cf. [16], Section 1.1.) 

For n=l, f(z) was analytic in the upper- and lower-half planes. In more than one 
dimension, /(z) need not be analytic, even though, for brevity, we_ still write it as a 
function of z rather than z and its complex conjugate z*. However, /(z) does in general 
possess a partial analyticity which reduces to the above when n = 1. Consider the partial 
derivative of /(z) with respect to z\ = Xk — iyu-, defined by 

2S t /=2|i = (55) 
dz* k dx k dy k 

Then / is analytic at z if and only if 3k f = for all k. But using our definition of /(z), 
we find that 
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2 * /(z) =^/_>l^ (x+Ty) - < 56 > 



It follows that the complex <9-derivative in the direction of y vanishes, i.e. 

/OO Qj 
dr (x + ry) 



»oo 



(57) 

if / decays for large |x| (for example, if / is a test function, as we have assumed). Equiv- 
alently, using 

2d k [6>(p • y) e 2 ^ z ] = 2d k [6>(p • y)] e 2 ^ z 

= . dfl(p-y) e 2^ P . z 

(58) 

= ip fc <5(p • y) e 2 ^ z 
= »p fc S(p • y) e 2 -P' x , 

we have for y ^ 

2y ■ 0/(z) = i / eTp (p • y) 5(p • y) e 2 ™ p ' x /(p) = 0. (59) 
Thus /(z) is analytic in the direction y. In the one-dimensional case, this reduces to 



dz* 



= Vy^O, (60) 



which states that f(z) is analytic in the upper- and lower-half planes. In one dimension, 
there are only two imaginary directions, whereas in n dimensions, every y ^ defines an 
imaginary direction. 

The multivariate AST is related to the Hilbert transform in the direction y (cf. [26], p. 
49), defined as 

1 f°° du 

H y f(x) = -PV / - /(x - uy), x, y G R n , y ^ 0. (61) 

(Usually, it is assumed that y is a unit vector; we do not make this assumption.) Namely, 
an argument similar to the above shows that 



/(x) ± z# y /(x) = 2 lim /(x ± iey), (62) 

ej.0 



hence 
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#y/( x ) = i lim [/(x - iey) - /(x + iey)]. 



(63) 



For n = 1 and y > 0, this reduces to the previous relation with the ordinary Hilbert 
transform. 

As in the one-dimensional case, f (x) is the boundary- value of f (z) in the sense that 

/(x) = lim [/(x + zey) + /(x - iey)]. (64) 
For real-valued /, these equations become 

/(x) = 2hm3?/(x + iey) 

~ (65) 
# y /(x) = 21im9/(x + iey). 

ej.0 

Two unresolved yet fundamental questions are: 

(a) For what classes of 'functions' (possibly distributions) can the AST be defined, apart from 
5(R n ); i.e., what is the domain of the AST? 

(b) Given a vector space H of 'functions' on R n for which the AST is defined, what is the 
range of the AST on HP That is, given a function F on C n , how can we tell whether F is 
the transform of some / G HI 

A necessary, though probably not sufficient, condition for F = / is that F satisfy the 
directional Cauchy-Riemann equation y • 9F(x,y) = 0. Complete answers to the above 
questions can be given in some specific cases: When / is a solution of the Klein-Gordon 
equation, then / must satisfy a certain consistency condition (the reproducing property of 
the associated wavelets). This condition, when satisfied by F, also guarantees that F = f 
for some / (cf. [16], Chapters 1 and 4). A similar result will be obtained in Section 4 for 
solutions of the wave equation in two space-time dimensions. The comments below apply 
to the general case and are, consequently, informal. 

The most direct way to find if F is the AST of some / is to construct f from F and then 
check that f = F. The first part has already been done formally, since / has been shown 
to be the boundary-value of /. Here we suggest an alternative method which can be used 
to construct f(p) instead of /(x). Assume that the Fourier transform is defined on H. If 
F(x, y) = /(x + iy) for some f (x) G H, then the 2n-dimensional Fourier transform of F 
is seen (formally) to be 

F(p,q)= f dxdye- 27 "( p - x+q ' y )F(x,y) 

= /(p) / dy6(p-y)e- 2 ^-^ ( 66 ) 

= /(p)<5(q-*p), 

where 5 is the AST, in Fourier space, of the Dirac measure 5(q). (This suggests that 
the AST, like the Fourier transform, exhibits some symmetry between space and Fourier 
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space.) 5 can be shown to be invariant under real rotations (q — ip — > P(q — zp), with 
R G SO(n)) and to transform under dilations as 

5(A(q-zp)) = A-^(q-zp), A^O. (67) 

Given p ^ 0, let R be a rotation such that p = |p| Re, where e = (1, 0, • • • , 0), and let 
k = |p| _1 P _1 q, so that q — ip = |p| R (k — ie). Then 



d(q-ip) = \p\- n d(k-ie) 



IPI 



/ dy^(yi) 

/ d yi e- 27T(1+lkl)yi [ dy 2 ■■■dy n e- 27Tt(k2y2+ - +knyn) 
Jo Jr.™- 1 



_ ^(fc2,---,fcn) 

~~ 27r|p| n (l + iki)' 

Let P : R n — > R and Q : R n — > R n_1 denote the projections Pk = fci and Qk = 
(k2, • • • , k n ). Then the numerator in the last expression is 

5(Qk) = SiQlpl^R- 1 ^)] = Ip)"" 1 SiQR^q), (69) 

hence 



^(q-»p) = o-,T ■ -oo-x v p^°- ( 7 °) 



2tt(|p| +iPR- 1 q) 

Together with eq. (66), this gives an explicit formal condition for F to be the AST of / 
and, if so, to determine /(p). When n = 1, 8 takes the simple form 

= afgrb- (71) 

3.3. Some Applications of the AST 

The AST is an example of a Windowed X-Ray transform, introduced in Section 2.1, with 
the window function 

h(r)* = —3 t-. (72) 

v ; 2tu(t-z) v ' 

(This window function is not "admissible" in the sense of Section 2.2, hence the recon- 
struction formula developed there fails. However, that reconstruction was based on some 
assumptions which may not be appropriate in general; see the comments below Eq. (88).) 
We now give two examples of the usefulness of the AST. An extensive use of this transform 
will also be made in Section 4. 
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Example 1 : Hardy spaces 

Suppose that /(p) vanishes outside of some closed convex cone V+. The cone V\_ dual to 
V + is defined as 

Vl = {ye R n such that p • y > Vp G V+}, (73) 

and it is clearly an open convex cone. Note that for y 6 7{, 6>(p • y) = 1 on the support 
of / (except at p = 0, which has measure 0), hence if / G L 2 (R n ) and y G V+, then 

/(x + iy)=/ rfpe 2 "P-^/(p) (74) 

and it follows (Stein and Weiss [27]) that /(z) is analytic in the tube domain 

T+ = {x + zy G C n such that y G (75) 

The set H 2 = {/ such that / G L 2 (V+)} is known as a Hardy space. Note also that / 
vanishes in the tube 



T" = {x + iy G C n such that -ye V^}, (76) 
since there p • y < for all ^ p G V+. Eq. (74) gives 

/(x) = lim/(x + zey), y G V[ , (77) 

ej.0 

which states that / is a boundary-value of /. Since /(z) is analytic, it may be regarded 
as a regularization of /(x) (the latter, being merely square-integrable, is a distribution). 
Eq. (77) can be viewed as a "reconstruction" of / from /, albeit a somewhat trivial one. 

Example 2: The Klein-Gordon Equation 

An important application of the AST is to signals that satisfy some partial differential 
equations. (In fact, it was in this context that the transform originated.) Suppose that / 
is a solution of the Klein-Gordon equation in R n , 

□/ + m 2 c 4 / = 0, (78) 

where 

D= d^_ c 2_d^_ _ c 2_^_ (79) 

is the D'Alembertian or wave operator for waves with propagation speed c. Here R n is 
interpreted as space-time, with x\ the time coordinate and (x2, • • • , x n ) the space coordi- 
nates, and m > is a mass parameter. This equation describes free relativistic particles 
in quantum mechanics. The limit m — > gives the wave equation, which will be discussed 
below. Define the solid light cone in Fourier space by 

V = {p G R n such that p 2 = pj - c 2 p 2 2 c 2 p 2 n > 0}. (80) 
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(Note that we are now using a Lorentz metric!) V is the union of the forward and backward 
light cones V+ and V_, where p\ > and pi < 0, respectively. Note that V± are convex 
but V is not. The fact that / satisfies the Klein-Gordon equation means that its Fourier 
transform /(p) is supported on the double mass hyperboloid 

fi m = {p G R n such that p 2 = m 2 c 4 } = 0+ U Q" , (81) 

where fi^j C V±. Thus / = /+ + / _ , where J 1 * 1 are distributions supported on f2±. Since 
C V±, an argument similar to that used for Hardy spaces shows that the corresponding 
solutions / ± (x) have AST's / ± (z) which are analytic in T ± and vanish in T T , where 

T ± = {x + zy G C n such that y G V±} (82) 
and V± are the cones dual to V±, which can be seen to be 

Vi = {y G R n such that y 2 = c 2 y\ - y\ y 2 n > 0, ± Vl > 0}. (83) 

Note that while V is a cone in Fourier space (i.e., p% is a frequency and p2, ■ • -p n are wave 
numbers per unit length), V' = V+ U VZ. is a cone in space-time. Technically, these two 
spaces are dual and should not be identified with one another. 

The AST of /, given by /(z) = f + (z) + /~(z), is therefore analytic in the double tube 
T = T + U T~, with T + and T~ containing only the positive- and negative-frequency 
parts of /, respectively. This "polarization" of frequencies is important because it makes 
it possible to reconstruct the solution / from / without approaching the singular boundary 
R n (y — *■ 0)- Eq. (74) shows that have the form 

;±(z)=/ dpet(pTa ± (p), zgT±, (84) 

where 

_ dp2 dp 3 --- dp n _ dp 2 dp 3 --- dp n 
P ~ " 2|pi| ~ 2 Cv /m 2 c 2 + p'j + ■ ■ ■ + p 2 n 

is the induced measure on Vt m and 

e z ± (p)*=e 2 -P' z , zgT±, P GO±. (86) 
The corresponding expression in the space-time domain, defined by 

e±(x')* = / dp e 2T*p-(— *'), (87) 

is a solution of the Klein-Gordon equation which can be shown to be a coherent wave- 
packet whose parameters z = x + iy have a direct geometric interpretation: x is a 
point in space-time about which is "focused" (i.e., e^(x') converges toward the point 
(x2, • • • , x n ) in space for times x[ < x\ and diverges away from it for times x[ > x\), and 
y is a set of homogeneous coordinates for the average velocity at which is traveling. 
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Furthermore, the invariant A > defined by A 2 = y 2 can be interpreted as a scale param- 
eter which, roughly speaking, measures the spread (resolution) of the wave packet at the 
instant of its maximal focus (x^ = x\) in its rest-frame (the coordinate system in which 
Di = 2/3 = • • • = Vn = 0). For small A, e^(x') is localized near x' k = Xk (k = 2, 3, • • • , n) at 
time x[ = xi, whereas for large A, it is spread out in space. (In Fourier space, on the other 
hand, it is A -1 which measures the spread.) The positive- and negative-frequency packets 
e+ and e~ are interpreted in quantum theory as particles and antiparticles, respectively. 
(This agrees with the usual observation that antiparticles "go backward in time." Cf. [16], 
Chapter 5.) 

Since the window function h(r) used in the AST has Fourier transform /i(£) = 6>(£) e -271 "^, 
it is not admissible in the general sense developed in Section 2.2; i.e., we have 



j^\h(0\ 2 = oo. (88) 

However, the rules of the game have changed. Eq. (88) was associated with a reconstruction 
formula which represents / as an integral of generalized wavelets parametrized by all of 
R n x R™, i.e., all x + zy with y ^ 0. This was acceptable when considering general 
functions /(x) in L 2 (R n ), since then we could define a representation of the affine group 
on such functions. But now we are dealing with a Hilbert space 7i of solutions of the 
Klein-Gordon equation, 

/(x)= ! dpe 2 ^a(p), (89) 



with the Sobolev norm 



dp\a(p)\ 2 (90) 



rather the L 2 norm used in Section 2.2. General affine transformations no longer map 
solutions to solutions, i.e. they no longer define operators on H. The mass m spoils the 
invariance of the equation under dilations. Only the subgroup V of translations together 
with Lorentz transformations (i.e., linear maps y i— > Ay which preserve the Lorentz norm 
y 2 ) maps solutions to solutions. V is called the inhomogeneous Lorentz or Poincare group. 
Recall that the measure used in the reconstruction formula of Section 2.2 was chosen to 
be invariant under dilations and rotations. Since dilations no longer define operators on 
H, this measure is no longer appropriate. Rather, we now expect to reconstruct / by 
integrating in y over the double hyperboloid 

fi A = { y g V such that y 2 = A 2 } = 0+ U (91) 

for an arbitrary fixed A > 0. Furthermore, we do not expect to integrate over all x e R n , 
since a solution is determined by its data on any Cauchy surface S C R n . For simplicity, 
take S to be the hyperplane x\=t for fixed t G R, though any Cauchy surface (spacelike 
(n — 1) -dimensional submanifold of R n ) will do ([16], Section 4.5). Thus consider the 
(2n — 2) -dimensional submanifold 
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a = {x + zy G T such that X\ = t, y G = a+ U cr_, (92) 

where y G in o±. a parametrizes all possible locations and velocities of a classical 
particle at the fixed time t, i.e. it is a phase space. A reconstruction formula has been 
obtained in the form 

/(x') = / ^(z)e z (x / )/(z), (93) 
where e z = on a±, a± is parametrized by (x 2 , ■ • -x n , y 2 , ■ • -y n )j an d 

d/i(z) = A(X, m)~ 1 dx 2 ■ ■ ■ dx n dy 2 ■ ■ ■ dy n . (94) 

A(X, m) is a certain constant related to the admissibility of the wavelets e z with respect 
to the measure dx 2 ■ ■ ■ dy n . Note that this differs from the usual construction of a solution 
from its initial data, which uses the values of both / and df jdx\ on S. The intuitive 
explanation is that the dependence of /(x + iy) on y G for fixed x G S, gives the 
equivalent "velocity" information. The independence of the reconstruction from the choice 
of Cauchy surface is due to a conservation law satisfied by solutions (cf. [16]). 

The above reconstruction formula bears a close resemblance to the standard representation 
of a function in terms of wavelets, for the following reason: In the hyperbolic geometry of 
spacetime, a moving object undergoes a Lorentz contraction, i.e. it shrinks in its direction 
of motion. Since y represents a velocity, eq. (93) expresses / as a linear combination of 
"wavelets" centered about all possible points in space (at the given time t) and in various 
states of compression. However, the analogy is incomplete since the e z 's can only contract 
and not dilate. (That is, they have a minimum width in their rest frames, determined by 
the choice of A.) Their contraction is due to Lorentz transformations rather than ordinary 
dilations of the form x 1— > ax, a 7^ 0. As noted earlier, the Klein-Gordon equation is not 
invariant under such dilations, due to the presence of m > 0. On the other hand, the wave 
equation (m — > 0) is invariant under dilations, hence the analogy with wavelets can be 
expected to be closer. This is the subject of our next section. 

4. Wavelets and The Wave Equation 

4.1. Introduction 

As explained in Section 1, a representation of solutions of the wave equation similar to 
that given for the Klein-Gordon equation in (93) should be of some interest in the analysis 
of naturally occurring signals, since it would automatically display much of their infor- 
mational contents. Unfortunately, the reconstruction formula in (93) fails when m — > 
because A(X, m) diverges in that limit. (The wavelet representation is no longer square- 
integrable in that limit.) This is probably related to the fact, well-known in quantum me- 
chanics, that the Klein-Gordon equation with m > has a very different group-theoretical 
structure from the wave equation. The symmetry group of the Klein-Gordon equation is 
the Poincare group V, while the symmetry group of the wave equation is the conformal 
group C, which contains the Poincare group as well as dilations and uniform accelerations. 
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The fundamental difference between massive particles (such as electrons) and massless 
particles (such as photons) is that the former can be at rest while the latter necessarily 
travel at the speed of light. It may well be that once the confer mal group is taken into 
account, an appropriate reconstruction formula can be found. In this section we confirm 
this hypothesis in two-dimensional space-time. 

4.2. Symmetries of the Wave and Dirac Equations in R 2 

In this subsection, we examine some group-theoretical aspects of solutions of the wave 
equation in R 2 . To simplify the notation, we choose the units of length and time such that 
the propagation speed c = 1. The wave equation then reads 



o = -n/(x,t) = (a 2 -a 2 ) f( x ,t) 

= (d x + dt)(d x -dt)f(x,t), 



(95) 



and we consider solutions which are possibly complex- valued. In terms of the light-cone 
coordinates 

u = x + t, v = x — t, (96) 

the equation becomes d u d v f = 0, hence the general solution has the form first given by 
D'Alembert, 

f(x,t) = f + (u) + f.(v). (97) 

f + (t + x) is a left-moving wave since it is constant on the characteristics x = xq — t, and 
similarly /_(t — x) is a right-moving wave. Later we shall find an appropriate family of 
Hilbert spaces 7i s (s > 0) to which f± will be required to belong, so we now write f± G 7i s 
without being specific. Note that we can let f± — > f± ± c, where c is a constant, without 
affecting /. This ambiguity will not be a problem since H s contains no non-zero constants. 
Hence the solutions are in one-to-one correspondence with the elements of the orthogonal 
sum V s = H s © 7i s , whose elements can be written in the vector form 

V>=(£)- (98) 
The wave equation can be restated as d v f + = d u f- = 0, or 

( dt -d x 9l J 4 )^(7A + 7A)^W = 0, (99) 

where 

Eq. (99) is known in quantum mechanics as the (two-dimensional, massless) Dirac equa- 
tion; i/j is a spinor, and 7t,7 x are Dirac matrices. The Dirac equation may be viewed 
as a particular system of first-order equations equivalent to the wave equation which is, 
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moreover, especially suited to the symmetries of the latter. The Dirac operator ^ is a 
"square root" of the wave operator in the sense that ^ 2 = — where I is the 2x2 iden- 
tity matrix. (This is related to the Clifford algebra associated with the Lorentz metric.) 
In more than one space dimension, there is no such simple relation between scalar-valued 
solutions / and spinor- valued solutions ip. 

We denote by V s+ and V s - the subspaces of V s with vanishing second and first compo- 
nents, respectively. Thus V s ± H s . 

A symmetry of the wave equation is a transformation which maps solutions to solutions. 
As symmetries can be composed and inverted, they form a group. We shall be particularly 
interested in "geometric" symmetries, which are induced from mappings on the underlying 
space-time that respect the wave equation. Some obvious ones are: 

Translations: For each (xo,to) £ R 2 , the map (x,t) — > (x + xo,t + to) induces a symmetry 
transformation f(x, t) — > f(x — Xo, t — to). This maps the right- and left-moving waves by 

T(u ,v ): f + (u)^f+(u-u ), f_(v)^f_(v-v ), (101) 

where u$ = xq +to and vq = — to- Hence translations can be made to act independently 
on the left and right parts of solutions. 

Lorentz transformations: For any real 6>, the map 

x — > £cosh6> + tsinh6>, t — > x sinh 9 + t cosh 9 (102) 

preserves the Lorentz metric x 2 — t 2 = uv. (It is a "rotation" by the imaginary angle 
—i9 and represents the space-time coordinates as measured by an observer moving with 
velocity — tanh6>.) This map has a particularly simple form in terms of the light-cone 
coordinates: 

u^e e u = Xu, f^e"% = A"V (103) 

Since this leaves the wave operator □ = —4d u d v invariant, it induces a symmetry on 
solutions. The simplest such transformation acts on f± by f+(u) — > f + (X~ 1 u), f-(v) — > 
f-(Xv). We shall need a more general induced transformation, given by 

L(A) : /+(«) - S+(X) MX-'u), /_(«) - S_(X) f_(\v), (104) 
where S± : R + — > C* = C\{0} (called "multipliers") must satisfy 

^(A" 1 ) = ^(A)" 1 , S±(XX') = S±(X) S±(X') (105) 

in order that Lorentz transformations form a group. (This means that S± are group 
homomorphisms.) Continuity in A > then implies that S+(X) = A - - 7 and S-(X) = A - - 7 
for some j, j' G C. 

The three-dimensional symmetry group of all maps T(uq, vq) L(X) is the restricted Poincare 
group Vq in 1 + 1 space-time dimensions. Note that Vq leaves invariant the subspaces V s ± 
of right- and left-moving waves. Since we shall be interested in irreducible representations 
of the symmetry group which characterize the complete wave equation, it is desirable to 
include a symmetry which mixes these two subspaces. 
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Space reflection: The map (x,t) — > (—x,t) is a discrete symmetry of the wave equation, 
corresponding to (u, v) (—v, —u). We take the induced mapping on solutions to be 

P: /+(«)->/_(-«), /_(„)-/+(-„). (106) 

Thus, P interchanges right and left waves, as desired. The idea that right and left be on 
equal footing is expressed as 

L(X)P = PL(A" 1 ), (107) 

which implies that 

S-(\- 1 ) = S + (\), (108) 

hence j' = —j. We call j the Lorentz weight of the solution. (In more than one space 
dimension, it is related to spin.) 

The group obtained from Vo by adjoining the space reflection is called the orthochronous 
Poincare group in the physics literature, since it still leaves the direction of time invariant. 
The full Poincare group V is obtained by further adjoining time reversal. However, the 
latter must be antilinear for reasons which need not concern us here (cf. Streater and 
Wightman [28]). The fact that leaves the direction of time invariant implies that the 
subspaces of positive-and negative-frequency solutions are invariant under it. To mix 
them, we introduce the following. 

Total reflection: The map (x,t) — > (—x,—t) is another discrete symmetry of the wave 
equation, corresponding to (u, v) — > (—u, —v). We take the induced mapping on solutions 
to be 

R: /+(«)->/+(-«), /_(«)->/_(-«). (109) 

Note that unlike translations, Lorentz transformations do not act independently on the 
left- and right-moving waves. This will be remedied by including the next symmetry. 

Dilations: Since the wave equation is homogeneous in x and t, it is invariant under the 
map (x, t) — > (ax, at), for any a^O. Equivalently, (u, v) — > (au, av). It suffices to confine 
our attention to a > 0, since dilations with a < can be obtained by combining D(a) 
with R. As with Lorentz transformations, we shall allow a multiplier M : R + — > C* in 
the induced mapping. Thus f(x,t) — > M(a) f(a~ l x,a~ 1 i), or 

D(a) : f+(u) - M(a) f+ia^u), f_(v) - M(a) /-(a^v). (110) 

In order that dilations form a group, we must have M(ct _1 ) = M(ct) -1 and M(aa') = 
M(a)M(a'). Again, continuity implies that M(a) = a~ K for some k E C, called the 
conformal weight of the solution. 

The following simple argument should convince the reader of the need to include non-trivial 
multipliers in eq. (110), i.e., to consider conformal weights other than zero: Suppose that 
f(x,t) is a solution with conformal weight k and let g(x,t) = (ad x + bdt) f(x,t), where a 
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and b are constants. Then g is also a solution of the wave equation, and it is easily seen to 
have conformal weight k + 1. In this way we can shift the conformal weight of a solution 
by any positive integer n by applying a homogeneous partial differential operator of order 
n with constant coefficients. (In three space dimensions, the electromagnetic potentials 
satisfy the wave equation in free space; the electromagnetic fields are obtained from them 
by applying first-order partial differential operators, hence are solutions but with higher 
weight.) 

Note that dilations commute with Lorentz transformations. Composing L(X) and -D(a), 
we obtain 

/+(«) - a~ K X~ j f + (X~ 1 a~ 1 u), /_(«) - a~ K X j f-iXa^v). (Ill) 

In order for the combination of dilations and Lorentz transformation to act independently 
on / + and /_, we must therefore require that j = k. We shall refer to k = j = —j' simply 
as the weight of /. Setting (3 = aX and 7 = a/A, we then obtain 

/+(«) - P~ K MP-'U), f-(v) -> 7" K /-(7"M- (112) 
Consequently, the semi-direct product R + x P acts on / by 

/+(«) - r K /+(/? _1 (u - uo)), /-(«) -> 7" K /-(T" 1 ^ - «o))- (113) 

This means that the group Go = R + x Vo generated by the continuous transformations 
T, L and D can be represented as a direct product Ax A of two copies of the affine group 
A acting independently on the right-moving and left-moving waves. We denote by Q\ the 
group obtained from Q$ by adjoining R, and by Q2 the group obtained by further adjoining 
P. Qi contains negative as well as positive dilations, but the signs of the dilations in the 
two components are equal. (They can be made independent by adjoining yet another 
discrete symmetry, namely (x,t) — > (t,x), but we resist the temptation.) Note that the 
subspaces of right- and left-moving waves are invariant under Q\ but not under £2- 
Since the affine group is closely related to wavelets, the decomposition Go A x A suggests 
that we apply separate wavelet analyses to /+ and /_ . Actually, we shall see that more can 
be done, due to the fact that the wave equation has another, quite unexpected, symmetry. 

4.3. Hilbert Structures on Solutions 

So far we have dealt exclusively with solutions in the space-time domain, where symmetries 
have a direct and intuitive meaning. We now wish to introduce inner products on solutions 
which make the above symmetry transformations unitary. For this purpose, it is necessary 
to venture into the Fourier domain and, later, also into the domain of complex space-time, 
invoking the Analytic-Signal transform introduced earlier. 

Let f(x,t) be a solution transforming under Q 2 with weight k. Representing f± formally 
by Fourier integrals, we obtain 

/oo /»oo 
dpe 2m P uf +{p)+ / dpe **ipv f_(p). (114) 
-00 J — 00 
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The symmetry operations defined earlier can now be represented in Fourier space. Trans- 
lations act by 

T(u ,v ) : hip) - e" 2 ^ f+ip), hip) - e~ 2 ^ /_(p), (115) 
Lorentz transformations act by 

L(A) : f+ip) - A^V+CAp), /_(p) - f-{\- x p), (116) 
space reflection acts by 

P : hip) - hi-p), hip) - Ui-P), (H7) 

total reflection acts by 

R - /±(p)->/±(-p) (118) 

and dilations act by 

D(a): f±(p)^a 1 -"f±(ap). (119) 

Let s = 23?(k) — 1 > and define 7i s to be the the Hilbert space of all 'functions' g(x) 
such that 

/oo 
rfpb|- s |^)| 2 <oo. (120) 
-oo 

In the terminology of Battle [3], this is the "massless Sobolev space of degree — s/2." The 
norm of a solution fix,t) (which may be identified with the corresponding spinor ip e T> s ) 
is then given by 

Ill/Ill 2 = IIIV'III 2 =II/ + H 2 +II/-I| 2 . (121) 

From the above actions it easily follows that the transformations T, L, R, P and D act 
unitarily on T> s , thus giving a unitary representation of Qi- Denote by Tif and 7i~ the 
subspaces of 7i s with support on [0, oo) and (— oo, 0], respectively. (Note that since s > 0, 
the value of /± at p = is unimportant; when s = 0, the origin has zero measure, and when 
s > 0, f± E H s implies /±(0) = 0.) The continuous symmetries T, L and D all preserve the 
sign of p, hence the group Q generated by them leaves invariant the subspaces of solutions 

p oo 

f+(x,t)= / dpe 2 ^ u f + ip) 
Jo 

f+M= [ dpe 2 ^ u f + ip) 

,Z (122) 

fl(x,t)= / dp e 2 ^ hip) 
Jo 

hix,t)= [ dpe 2 ^ v hip), 
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where the superscipts denote positive- and negative-frequency components, as in Section 
3. The complete solution is / = ft + fT + fZ + ft- Note that f-(p) with p > 

I i I 

represents a negative-frequency component since v = x — t. Thus ff e Hf but f_ eHf. 
The decomposition H s = Hf © Ti~ therefore gives a corresponding decomposition 

V s = Vf + © VJ + © Vf_ © V~_ , (123) 

where Vf + = Hf © {0} and Vf + = {0} © Hf. Let £>± = X>± + © Vf_. The restriction of 
the representation of Q 2 to Go leaves invariant all four subspaces V s T a (a, r = ±), hence 
it is reducible. When R is included, only the subspaces T> s+ and V s _ remain invariant. 
When P is further included, no invariant subspaces remain and we have an irreducible 
unitary representation of Q2 on T> s . This shows that the discrete symmetries P and R 
serve to 'weave' the four representations on T>f±, T>f ± (which are associated with the 
affine group rather than the wave equation) into a single representation characterising the 
wave equation as a whole. This was the purpose for which they were introduced. 
However, distinct values of k are not significantly different at this point. For if / has 
weight k, let k' G C with s' = 2^R.(k') — 1 and 

9±(p) = \p\ K '- K f±(p)- (124) 

Then the corresponding solution g(x, t) has weight k' and, furthermore, the map / — > g is 
unitary from T> s onto T> s ' ■ This shows that the UIR's with any two values of k are unitarily 
equivalent. If no more could be said, we might as well restrict ourselves to a single value of 
since solutions of arbitrary weight can be obtained by the above method. (Recall that 
we could shift the weight up by a positive integer by applying a differential operator; the 
above generalizes this process to arbitrary shifts.) Actually, it turns out that for certain 
special values of k, the group of symmetries can be significantly enlarged, and distinct 
values of k then give unitarily inequivalent representations of the larger group. Unlike A, 
however, the larger group no longer acts simply on the Fourier space; instead, its natural 
domain of action is the complexified space-time associated with the AST. For this reason 
it becomes necessary to re-express the norm of V s in terms of the AST / of /, as will be 
done next. 

4.4. Norms in terms of Analytic Signals 

We wish to give the norms defined above expressions which are local in space, that is, 
involve only values of f(x,t) at an arbitrary fixed time t, say t = 0. Since a knowledge of 
both f(x,0) and d t f(x,0) is required to determine /, we cannot expect such expressions 
to characterize the complete solution /. On the other hand, the positive- and negative- 
frequency parts f ± of / satisfy the first-order pseudo-differential equation 

-id t f ± = ±^T&if ± 1 (125) 

hence are determined by their initial values / ± (x,0). (Thus, instead of using f(x,0) 
and dtf(x,0) as initial data, we use the symmetric pair f ± (x 7 0).) We therefore consider 
separately the norms on Hf and H~ , given respectively by 
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f+|||2 _ 



/oo 
dp\p\~ 
-oo 
poo 

/ dp p- s 

JO 

X) 

dp \p\~ 



-oo 
oo 



dp p 



~o( P )\U(p)\ 2 + o(-p)\f-(p)\ 2 ] 
\U(p)\ 2 + \f-(-p)\ 

o(-p)\U(p)\ 2 + o(p)\f-(p)\ 2 
\h(-p)\ 2 + \f-(p)\ 2 



(126) 



The separation of positive and negative frequencies is characteristic of the AST; hence let 
us consider the AST of the complete solution /: 



1 f dT 

f(x + ix', t + it') EE — / —— f(x + TX', t + Tt') 
27TZ J — oo T I 

1 f°° d-r 
= — J_ — [/+(« + TU') +f-(v + TV')] 

= f+(u + iu') + f_(v + iv'), 



(127) 



where the imaginary parts of the space-time and light-cone coordinates are related by 
v/ = x' + t',v' = x'- 1' and /± are the AST's of f±. Recall from Section 3 that the AST's 
of f ± (x, t) are obtained by restricting / to the forward and backward tubes where 
±£' > \x'\. This suggests that the AST is well-suited for the analysis of the above norms, 
since it naturally breaks / up into its four components. 

The symmetry operations on real space-time extend to the complexified space-time by 
C-linearity, and these extensions induce transformations on / in the obvious way. For 
example, 



T(u ,v ): f+(u + iu') -> f + (u-u + iu'), f-(v + iv') -> f+(v - v + iv'). (128) 



(Note that the parameters of the symmetries are still real; e.g., we do not consider complex 
translations.) 

Since only the combinations (x + ix') ± (t + it') enter into /, it suffices to set x' = and 
t = 0. This is called the Euclidean region in quantum field theory, since the Lorentzian 
metric x 2 — t 2 becomes x 2 + t /2 , which is Euclidean. (See Glimm and Jaffe [10].) Note that 



/>, it') = f+(x + It') + f_( X - It') EE f+(z) + f-(z*) 



(129) 



is now presented as a sum of analytic and anti-analytic functions of z = x + it', as befits 
a solution of the "analytically continued" wave equation 



d 2 d 2 
+ 



f(x,it')=0, 



dx 2 ' dt' 2 

which states that f(x, it') is harmonic. In terms of the Fourier transforms, 
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(130) 



/oo 
dp 9(pt') e 2 ^ 2 f+(p) + e(-pt') e 2 ^ /_(p) | . 
-oo 

Hence the initial values of are given by 

F ± (x) = f±(x, 0) = lim f(x, it') = f(x, ±i0) 



(131) 



dp e 



±t'->0+ 
2-Kipx 



e(±p)f+(p) + e(TP)f-(p) 



(132) 



from which 



F ± (p) = 9{±p) U(p) + 6{tp) /_(p). 



(133) 



Returning to the norms, consider first the simplest case s = (i.e., k = | + ifx, with // 
real) : 



f±|||2 



(134) 



dp e(±p)\f + (p)\ 2 + 9(Tp)\f-(p)\ 2 

oo 

/oo 
-oo 

/oo 
-oo 

by Plancherel's theorem. We have therefore proved 

Theorem 4. The norms in the positive- and negative-frequency subspaces ofV can 
be be written as 

/oo 
efe|/(z,±z0)| 2 , 
-oo 

and the complete norm in T>q is therefore 



(135) 



dx 



|/>,z0)| 2 + |/>,-z0)|' 



(136) 



This is not quite local in space since it involves the boundary values 
f(x, ±i0). Rather, it is 'local' in the complex space generated by the AST. We refer to 
this property as pseudo-locality; it stands in the same relation to locality as the pseudo- 
differential equation (125) stands to the wave equation. Since f(x, 0) = f(x, iO) +f(x, —iO) 
and f(x, it') may be regarded as a regularized version of f(x, 0), the term 'pseudo-locality' 
seems appropriate. 

Next, fix k E C with s > 0. Can we obtain pseudo-local expressions for the norms with 
s > 0? We have 



27 



r±|||2 



= r dp\ P \- \e(± P )\h(p)\ 2 + 8(TP)\f-(p)\ 2 } 

J — oo L J 



-oo 

-OO 



poo 

= dpp- 
Jo 



\U(±p)\ 2 + \f-(TP)\ : 



(137) 



Theorem 5. For arbitrary s > 0, the norms \\\f ± \\\ s in Vf have the following pseudo-local 
expressions in terms of the restrictions of f to the Euclidean region: 



*±|||2 



oo /*oo 



oo JO 



= N S I dx I dt' t' s 1 \f(x,±it')\ 2 , 



where 



Hence the norm in T> s is given by 



N„ = 



F(s) 



(138) 



(139) 



oo poo 

2 - N s I dx dt'\lf\ a - 1 \f(x,ilf)\ 2 

oo J — oo 



Proof: We prove the theorem for |||/ + |||s. The proof for 
f(x, it') can be written as an inverse Fourier transform 



hence by Plancherel's theorem, 



(140) 

is similar. By eq. (131), 
(141) 



/oo poo 
dx\f(x,zt')\ 2 = / dp eipt^e-^'iMptf + ei-pt^e^'if-iprf 
-oo J —oo 



-oo 
"OO 



/oo 
dp e( P t') e-*^' \\h(p)\ 2 + \f-(-P)\ 2 • 
-oo L J 



(142) 



When t' > 0, then p > in the last integral. Integrating over t' with the factor t' s , 
exchanging the order of integration on the right-hand side and using 



dt' t 



i -/V s - 1 -Airpt' _ 



= (4np)- s T(s) = Nr L p 



1 „S 



(143) 



for s > 0, we obtain 



/oo poo poo 

dx dt't' s - 1 \f(x,zt')\ 2 = dpp~ s \f + (p)\ 2 + \U-p)\< 
-oo JO JO 1 



(144) 



?+|||2 



as claimed. 
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4.5. The Wavelets e z and their Mother 

Here we show that Theorem 5 establishes a resolution of the identity I s in V s (s > 0) in 
terms of wavelets e z G T> s parametrized by z = x + it' ', t' ^ 0, in the manner described 
in Sections 2.2 and 3.3. We shall define e z as an element of T> s whose inner product with 
/ G T> s is the value of the AST / at (x, it'), i.e., 

f(x, it') = (( e z , f)) s = ( e z+ , f+) a + {e z _,f-) a , (145) 

where ( •, • ) s and (( ■, -)) s denote the inner products in H s and V s , respectively. Since e z 
is to belong to V s , it must itself be a solution. Eq. (131) shows that its components in 
Fourier space are 

e z + (p) = 0(pt') \p\ 8 e- 2 ^* 

e z .(p) = e(-pt') \p\ s e~ 27Tipz . (146) 



(Recall that according to our convention, (g, f) s is anti-linear in g.) Note that e z _ = e z * + 
as elements of TL S , and that e 2 _ and e z+ are orthogonal not only as elements of the direct 
sum T> s = T> s+ © T> s _ (since e 2± G T> s ±) but also as elements of H s . General solutions 
clearly do not share this property. The positive- and negative-frequency components ef 
of e z are obtained by choosing t' > and t' < 0, respectively: 



Hence 



For t' > 0, we have 



e+ if t' > 
e~ if£'<0. 



e+ + (p) = 6{p) p s e~ 2lTipz * , t' > 
e+_(p) = 6»(-p)|p| s e- 2 ^ z , t'>0 
e 2 - + (p) = ^(-p)H s e- 2 ^, f<0 
e;_(p)=%)p s e- 27 ^ 2 , f<0. 



I +|||2 _. II + ||2 I II + ||2 
\ e z Wis ~~ ll e z+lls ~r ll e z— lis 



/oo 
-oo 



Since e 2 _ = e z *+ in 7Y S , we have for t' < 
hence 



(147) 



(148) 



9(p)e-^' +9(-p)e^'] {ug) 

' + !) 
(4^) s+1 ' 



, ll9 2r(s + i) , s 

ej ;= , LnJi . ( 15 °) 
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lll«.lll! = IWIi; + IK-|||J = ^j5^ wf/o, (i5i) 

since one of the vectors ef vanishes. This proves that e z £ V s as claimed. 

Theorem 5 can be given a suggestive form in terms of the orthogonal projection P(z) onto 
the one-dimensional subspace of e z in V s . For t' 7^ 0, 

|/0r, it')\ 2 = I (( e 2 , /» s | 2 = \\\e z \\\ 2 S « /, P(z)f)) s = (( /, P(*)/». • (152) 

Substituting this into eq. (140) and simplifying, we obtain 

^ Foe ^ Foe % ^ ^ P(Z)f))s = ^ ^ /))S ' (153) 

Upon polarization, this proves the weak operator identity 

poo poo 

*L*'Li* p(z)=l -' (1B4) 

giving a resolution of the identity in T> 8 in terms of the projections P(z). Note that 
dxdt'/t' is just the left-invariant measure on the affine group. We shall see below that 
x and t' can indeed be used to parametrize the affine group, with t' < corresponding 
to negative dilations. The fact that eq. (154) involves the left-invariant measure on the 
affine group shows that the weight function It'] 5-1 in the original measure only served as a 
normalization factor for e z , so as to make e z+ and e z _ analytic in z* and z, respectively, 
which in turn makes ((e z+ ,f)) s = f+(z) and ((e z _,f)) s = f-(z*) analytic in z and z*, 
respectively. 

The wavelet e z can be regarded as a regularization of the point x in space, with \t'\ as a 
measure of its diffusion. Hence e z transforms naturally under those symmetry operations 
which leave the Euclidean region invariant. These consist of space translations T(x ), 
dilations D(a) and the reflections P and R. Clearly the Euclidean region is not invariant 
under Lorentz transformations and time translations. (It is not difficult to extend the 
family of wavelets to one parametrized by general points (x + ix 1 , t + it') in complex space- 
time, with x' 7^ rkt' . This extended family is then invariant under Q2, and the above inner 
products are obtained by simply integrating over the submanifold defined by x' = t = 0. 
Then the invariance of the inner product under the full group Q2, known to hold in the 
Fourier domain, can be made manifest in the complex space-time domain by introducing 
conserved currents and using Stokes' theorem. Cf. [16], Section 4.5 for a related treatment 
of the Klein-Gordon equation.) 

To see how e z transforms under a symmetry operation, we need only apply eq. (145). 
Thus for translations, 



(( T(x ) e z , f)) s = (( e 2 , T(-x )f)) s = (t(-x ) /) (x, it') 
= f(x + xq, it') = (( e z+Xo , f)) s , 



(155) 
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hence 



T(x Q ) e z = e z+X0 . (156) 

Similarly, 

((D(a) e z J)) s =((e Zl D(a- 1 )f)) s 

= a K f(ax, iat') = a K (( e az , f)) s , 

hence 



(157) 



D(a) e z = ct K e az , a > 0. (158) 
Since P* = P' 1 = P and (Pf)(x,it) = f(-x,it'), we have 

Pe z =e. z ». (159) 
Finally, R* = R~ l = R and (Rf)(x, it') = f(-x, -it 1 ) implies 

Re z = e_ z . (160) 

Thus T(xo) and D(a) generate a single copy of A (as opposed to Go ~ A x A), and P 
and R extend this to include space reflections and negative dilations. We may therefore 
begin with a single basic wavelet, say 4> = G T> s+ (left-moving wavelet with z = i). Its 
components in Fourier space are then 

4>+(p) = 6(p)p s e~ 27Tp , 4>_(p) = 0. (161) 

For z = x + it' with t' > we obtain all the left-moving wavelets by applying T, D and R 
to (p: 

z+ ^ \ J (162) 

e;. + = (0 K T(x)RD(t')(f>, 

and the right-moving wavelets are obtained in the same way from Pep = ej_. Note that 
the choice of <p depends only on s and not on the imaginary part of k. We may regard eqs. 
(162) as providing a construction of the UIR with weight k by giving all of its wavelets. 
Modified versions of the left component 0+ of (/> have already appeared in the literature, 
in connection with the representation theory of the affine group. Namely, if we absorb the 
Sobolev weight function |p| _s into the elements g of 7i s by defining G(p) = \p\~ s ^ 2 g(p) 
(so that G G £ 2 (R, dp)), then cf> becomes 

&+(p) = 9(p)p s/2 e- 27rp , &_(p) = 0. (163) 

The function $+(p) with s = 1 first appeared in the classical papers of Aslaksen and 
Klauder [1], where it was used as a 'fiducial vector' to obtain a 'continuous representation' 
of the affine group. Since those papers contain the first instance of what is now called 
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continuous wavelet analysis, we suggest that 4> with s — 1 deserves to be called the Mother 
of all Wavelets. For s = 1, 2, 3, • • • , &+(p) also appears in the work of Paul [24] in connection 
with representations of the affine group and their extension to SX(2,R) (see below). It 
must be noted, however, that in our case the wavelets are prescribed by the problem at 
hand (solving the wave equation and extending the solutions to complex space-time, via 
the AST) rather than chosen arbitrarily as a convenient family of functions to be used in 
expansions. This is further discussed in Section 5. 

Having established a resolution of unity in terms of the wavelets e z , let us now investigate 
the associated reconstruction (cf. Section 2.2). The formula corresponding to eq. (22) is 

f{x,t) = N s [ dx^t'^tW 3 - 1 e Xl+u[ (x,t) f(x ui t[). (164) 
Jr. 2 

This gives an expansion of an arbitrary solution / e T> s in terms of the wavelets e Zl , this 
time expressed in the real space-time domain. Let us therefore compute the left-moving 
wavelet e Zl + (x, t). For simplicity, choose x\ = and t\ > 0; the other cases (x\ ^ 0, t[ < 
and right-moving wavelets) can be obtained easily by using T(xi), R and P. By eq. (146), 

= f 3 dppV**<«+*i> = - ^ (S + \L 1 - (165) 
Jo (27r(r 1 — iu)) 



Hence 



\e lt > + (x,t)\ 2 = Eii+l£ _ (166) 

At time t, this solution is localized in space around x = —t, and its width is proportional 
to t^. As t'i — > 0, the wavelet becomes an infinitely sharp spike at x = —t. This is not 
surprising, since t[ acts as a scale parameter for the affine group. Similarly, e Xl+ i t ' i (x, t) is 
centered near x = x\ — t. 

A direct and important interpretation of t' is that 1/t' is proportional to the average 
frequency v s ±, or color, in the frequency spectrum of e z ±. This can be seen most easily 
by noting that the function p s exp(— 2npt') has a maximum at p = s/(2irt'). A more 
precise argument is based on the quantum mechanical notion of expectation values, where 
l^z±(p)| 2 is viewed as an (unnormalized) probability distribution for p with respect to the 
measure \p\~ s dp. Remembering that the frequency in D s± is ±p, this gives 

JZo\p\~ s dp (±p)\e z ±(p)\ 2 ^ f 1 d\. „ ||2 s + 1 



Us± s "7_°ih-^i4 ± (p)p = ± K-tom) log l|e - L = i^- (167) 

(Note that z/ s+ = z/ s _ = z/ s and, as expected, v a has the same sign as t'.) Similarly, one 
computes the standard deviation in the frequency to be 



a "- = i/(-S5f) log l' e -"-W- (168) 
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Note that eq. (154) assumes an especially simple form in terms of the variables (x, v s ): 

[ dx f dv s P[x + i^-] = I s , (169) 
8 + 1 J.oo J.oo V 47rz/ s y 

which gives a resolution of the identity in T> s in terms of the orthogonal projections onto 
the wavelet subspaces parametrized by initial location and color, with all locations and 
colors given equal a priori probability since the measure is Lebesgue! (The appearance of 
Lebesgue measure is not an accident. The variables (x, v s ) are phase-space coordinates, 
also famous in symplectic geometry as Darboux canonical coordinates, and dx dv s is the 
corresponding Liouville measure. For a similar analysis of the wavelets associated with the 
Klein-Gordon equation, see ref. [16], Section 4.4.) 

We are now in a position to answer a question which was posed in a more general context 
in Section 3.2: Given a function F(x,t'), how can we tell whether F = f(x,it') for some 
/ G T> s ? Suppose this were the case for some positive integer s. From the polarized version 
of eq. (140), i.e., 

/OO nOO 
dx / dt'^-^g^it')* f{x,it'), (170) 
-oo J — oo 

we obtain by choosing g = e Zl with z\ = X\ + it[, z = x + it' : 

/OO POO 
dx / dt' {t'r 1 (( e Zl , e z )) s f(x, it') . (171) 
-oo J — oo 

The function 



K(zi,z) = ((e Zl ,e 2 



(172) 



is a reproducing kernel for the function space T s = {/ such that / G T> s }. A computation 
similar to that for |||e z ||| s gives 



{2TT\t' 1 +t'\) 



s+1 



1 



. X\ — X 

l t[T¥ 



-s-l 



(173) 



The integral operator defined by K acts as the orthogonal projection from L 2 (|t'| s_1 dx dt') 
to JF S . The given function F is therefore the AST of some / G T> 8 if and only if is satisfies 
the consistency condition 



/oo poo 
dx / dt'\t'\ s - 1 K(z 1 ,z)F(x,t') 
-oo J — oo 

(cf. [16], Chapter 1 for an exposition of this general idea). 



(174) 



4.6. Extension to SL(2, R) 

We are at last ready to extend the symmetry group to a larger one that will 'break' the 
unitary equivalence of the UIR's of A with different weights. We do so by showing that 
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for 2k = 1,2,3, the unitary action of A on each of the subspaces V s ± extends to 
the group of real 2x2 matrices of unit determinant, and that this extension preserves 
the positive- and negative-frequency subspaces T>^± and T>~ ± and is therefore irreducible 
when restricted to these subspaces. The extension can be implemented most easily on 
the AST's of solutions in V s . We concentrate on V s+ ks H s for simplicity and denote its 
elements by g rather than /+ to avoid a proliferation of indices. 

The action of the affine group on 7i s is induced by the map u — > au+[3 on R, which is but a 
very special case of an order-preserving diffeomophism h : R — > R; i.e., an invertible map 
h such that both h and h~ x are C°° and dh/du > 0. Under composition, the set of all such 
maps forms a group Diff + (R), sometimes called the pseudo-conformal group, of which 
A is a subgroup. Note that whereas it takes only two parameters to specify an affine 
map, an infinite number of parameters are needed to specify a general diffeomorphism. 
Nevertheless, for k = \, the action of A on Hq = L 2 (R) can be extended to Diff + (R) by 



which is unitary and gives a representation since A{hi) A{Ji2) = A{h\ o /i 2 ). However, 
this representation is problematic since general elements of Diff + (R) do not preserve 
the positive- and negative-frequency subspaces of Hq. To see this, recall that g is the 
boundary-value of its AST g(z) and g = in if g G Hf. Hence elements of Hf are 
boundary-values of functions analytic in C^. When composed with a diffeomorphism h~ x 
which is not the boundary-value of a function analytic in C 1 * 1 , the result cannot be the 
boundary-value of a function analytic in C^. (Diff_|_(R) plays important roles in string 
theory and quantum field theory; there, the fact that it mixes positive and negative fre- 
quencies leads to a difficulty known as the "conformal anomaly") Since the decomposition 
into positive- and negative-frequency components is fundamental to our approach, the 
representation theory of Diff + (R) will not be pursued further. However, the above discus- 
sion suggests that we confine ourselves to diffeomorphisms of R which extend 'naturally' 
to C and whose restrictions to C ± are holomorphisms (analytic diffeomorphisms) of C 1 * 1 . 
That is, we consider symmetry operations on H s which are induced from geometric maps 
of the complex space C + U C~ associated with the AST rather than geometric maps of the 
real space R. Since only the analytic functions g(z) or their boundary-values g(x±i0) en- 
ter into the pseudo-local inner products of H s , we might even get away with maps which 
have singularities on R. Thus we look for the group of holomorphisms of C ± . Now it 
is well-known that the group of holomorphisms of the Riemann sphere C = CU{oo} is 
SL(2, C), the set of all complex 2x2 matrices with unit determinant, acting by fractional- 
linear (Mobius) transformations. The subgroup preserving C^, hence also the one-point 
compactification R = R U {oo}, then consists of all such reai matrices, and is denoted 
£X(2,R). We shall write G = 5X(2,R) throughout this subsection. An element of G, 
being a diffeomorphism of R rather than R, will have a singulariy in R (the inverse image 
of oo) unless it has oo as a fixed point. 

We begin by constructing two inequivalent UIR's of G which extend the UIR's of A on 
H^. The matrix 



defining 




(175) 
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a b 
c d 



ad-bc=l, (176) 



acts on C + U C = C\R by the fractional-linear transformation 



az + b _ 

a(z) = ; = w. 177 

' cz + d ' 



(This action will be derived below.) From the easily derived identity 



z — z 



w-w* = , (178) 

\cz + d\ 2 

it follows that a leaves C ± invariant. Note that the boundary map a(u) on R has a 
singularity at u = —d/c, in agreement with the above discussion. Since da(u)/du = 
(cu + d)~ 2 , the induced map on boundary-values in Hq is, according to eq. (175), 

A(a~ 1 ) : g(u) - \cu + d\~'g . (179) 

\cu + d J 

The singularity at u = —d/c poses no problem, since C°° functions of compact support 
g G T>(R) vanish at infinity and £>(R) is dense in Hq. The set of cr's whose boundary maps 
are non-singular is, in fact, just the affine group A, since c = implies 

cr(u) = a 2 u + ab = au + f3, (180) 

giving the relation between the parameters of A and those of G. To represent A(ij~ 1 ) 
on g(z), we first replace the multiplier \cu + <i| -1 with (cu + d) -1 , which also leads to a 
representation and, unlike the former, can be extended analytically. This does not affect 
the norm as z — > u G R. Thus define 

B(a~ 1 ) : ~g(z) - (cz + d)^~g {^j^j ■ (181) 

For z G C 1 * 1 , we obtain the actions on 7i^, which are unitary by the same argument as 
used for Diff + (R). This gives inequivalent UIR's of SX(2,R) on H^. Comparison with 
Bargmann's classification shows that these coincide with the two representations of the 
"mock discrete series" with s = (Lang [20], pp. 120 and 123). 

We can arrive at the action of G on C 1 * 1 (eq. (177)) by the method of cosets. Consider the 
subroups T and H of G given by 



(182) 



t (u) = ^ ^ such that u G r| 

H = \h(c,d) = (^^ ^ such that d^oj. 

Then T is isomorphic to R, and H is isomorphic to A. Moreover, any element of G with 
d j£ can be written uniquely in the form 
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-(c 5) = G 2)='W»M. (183) 

Hence we may coordinatize G (except for a set of Haar measure zero) by the product of 
sets TH. (This does not mean that G is a direct product since T and H, as subgroups of 
G, do not commute.) The space U = G/if of right cosets is then parametrized by u G R, 
as follows: If d ^ 0, then 

a# = {a/c such that /c G -ff} = {t (b/d) h(c, d) k such that k G H} = t (b/d) H = H b/d . 

(184) 

The set of matrices with d = forms a single coset, which will be denoted by H^. Thus 
U « R. Now G acts on £/" by left multiplication, and this translates into an action of G 
on R as follows: If u ^ oo, then 

(: s)*.=(: 9 (J :^)*=^ ^> 

with ct('u) as given by eq. (177). The action on is similarly found to be 

^ ^ 00=J ff a/c . (186) 
(If c = 0, then is invariant under a.) 

Thus we see that the action of G on R (more precisely, on R) can be derived from the action 
of G on itself. Similarly, the representation of G on £ 2 (R) obtained above can be viewed 
as a group-theoretical construction: The Lebesgue measure du on R can be extended to 
U = R by letting {oo} have measure zero; this extension du is quasi-invariant under the 
action of G, and functions in L 2 (R) = L 2 (U,du) may be interpreted as functions on G 
which are constant on each coset H u . The above representation can then be interpreted 
as being "induced" (Lang [20], chapter 3) from the identity representation of H. Since 
H is not compact, non-zero functions on G which are constant on the cosets H u are not 
square-integrable on G. This proves that the above representation is not square-integrable 
on G. That means that we cannot use it to obtain a resolution of the identity in L 2 (G). 

We now construct some other representations of G, with weights k other than ^, which 
do turn out to be square-integrable on G. Let k be real, so that s = 2k — 1. Define the 
action of G on H s by 

Bio- 1 ) : g(z) - (cz + d) — 1 ~g • (187) 

We must have s G Z in order to preserve analyticity. For k = |, this coincides with the 
earlier action. Note that for a = b = —{3/y/a, c = and d = y/a we obtain 

Bio- 1 ): g(z)^a- K g(a-\z-(3)), (188) 
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which shows that the action restricts to that of the affine group on H s - 

To show that B(a~ 1 ) acts unitarily, we need to prove that it preserves the pseudo-local 
norm introduced in Theorem 5. There we found that for G replaced by its subgoup 
A, we needed s > 0. Since s must now, in addition, be an integer, we therefore have 
2k — 1 = s = 1, 2, 3, • • •. We must show that 

/ dPz\Q(z)\ a - 1 \cz + d\- 2a - 2 g(^^\ = [ d 2 z\Z(z)\ s - 1 \g(z)\ 2 , (189) 
Jc \cz + dj J c 

where d 2 z denotes Lebesgue measure in C. Let w = <j(z) = (az + b)/(cz + d). Then 
dw/dz = (cz + d)~ 2 , hence 

d 2 w = \cz + d\~ 4 d 2 z, and S(w) = \cz + d\~ 2 S(z), (190) 

where the second equality follows from eq. (178). Thus 

/ d 2 z\Z(z)\ s - 1 \cz + d\- 28 - 2 \~g(w)\ 2 = [ d 2 w\Z(w)\ 8 - 1 \~g(w)\ 2 =\\g\\i (191) 
Jc Jc 

proving the result. The unitary representation on H s decomposes into a direct sum of 
UIR's on Tif. All these representations are inequivalent, and they are known collectively 
as the discrete series (Gelfand et al. [9], Lang [20]). 

The action of SL(2, R) on the wavelets e z+ G V s+ ks H s is easily computed. By the 
unitarity of B(a), we have 



(B(a)e z+ ,g) s = ( e z+ , B{a~ l ) g ) s = {B{cj- l )~g) (z) 

= (cz + d)-*- 1 ~g(a(z)) = (cz + rf)" 5 " 1 ( e a(z)+ , g" 

hence 



B(a) e z+ = (cz* + d)"*" 1 e a{z)+ . (193) 

The representations of the discrete series are square-integrable over G. The subgroup of 
G which leaves z = ±i invariant is SO(2), hence ~ G/SO(2). Since SO (2) is compact 
(unlike its counterpart H for s = 0), the norms in Ti.f can be rewritten as integrals over 
all of G rather than just the homogeneous spaces C 1 * 1 . 

Returning to the wave equation, we obtain mutually inequivalent UIR's of what we shall 
call the restricted conformal group Co = GxGon , V~ + , V~_ and (s = 0, 1, 2, • • ■). 
When the total reflection R is included, we obtain inequivalent UIR's of the resulting group 
C\ on T> s+ and T> s _. When the space reflection P is further included, we obtain a single 
set of mutually inequivalent UIR's of the resulting group C2 on V s . As they did for Qq, the 
reflections unify the four subspaces Vf ± into a single one representing the wave equation 
as a whole. 
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5. Concluding Remarks: Dedicated Wavelets 

The wavelet analysis developed in [16] for the Klein-Gordon equation has the interesting 
feature that the wavelets are "dedicated" to the equation rather than being merely a con- 
venient set of functions to be used in expansions. (This is somewhat reminiscent of the 
situation in the spectral theorem, where expansions are customized to a given operator.) 
The reward for such dedication is that symmetry operations (such as translations, rota- 
tions, Lorentz transformations, and even time evolution) take wavelets to wavelets. This 
has the practical consequence of making the description economical and precise. For ex- 
ample, wavelets obtained by taking tensor products of one-dimensional wavelets cannot be 
rotated; consequently, a function consisting of but a few wavelets in one coordinate system 
is represented (inefficiently) by a combination of many wavelets in a rotated coordinate 
system. Similar considerations apply to the dedicated wavelets associated with the wave 
equation in R 2 developed in the last section. But since there is now only one space dimen- 
sion, the results are somewhat less dramatic: There are only two directions, left and right, 
giving rise to the labeling e z+ and e z _. We believe that the results of Section 4 generalize 
to d > 1 space dimensions, where the set of directions is parametrized by S"^ -1 . (The case 
d = 1 is degenerate since S° = {±1} is disconnected.) Work on this is in progress. 
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